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We put forward a new, versatile and highly-scalable experimental setup for the realization of 
discrete two-dimensional quantum random walks with a single-qubit coin and tunable degree of 
decoherence. The proposed scheme makes use of a small number of simple optical components ar- 
ranged in a multi-path Mach-Zehnder-like configuration, where a weak coherent state is injected. 
Environmental effects (decoherence) are generated by a spatial light modulator, which introduces 
pure dephasing in the transverse spatial plane, perpendicular to the direction of propagation of the 
light beam. By controlling the characteristics of this dephasing, one can explore a great variety 
of scenarios of quantum random walks: pure quantum evolution (ballistic spread), fast fluctuat- 
ing environment leading to a diffusive classical random walk, and static disorder resulting in the 
observation of Anderson localization. 

PACS numbers: 03.67.-a,03.67.Ac 03.65.Yz 
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I. INTRODUCTION 

The underlying principles of many quantum informa- 
tion protocols can be traced back to the concept of quan- 
tum random walk (QRW), which since its first descrip- 
tion have become a fundamental paradigm in quantum 
science [l|,[2j]. The idea of QRWs was originally conceived 
by Aharonov et al. [3] as an extension of the well-known 
classical random walk (CRW) [4]. The main distinguish- 
ing feature of a QRW, compare to a CRW, is the possibil- 
ity of interference between the multiple paths that can be 
simultaneously traversed by a quantum walker, enabling 
thus a faster spreading of the uncertainty of location of 
the walker than in the classical case [5|, [6| . 

The temporal evolution of a quantum system, such as 
a QRW, depends on the presence, and specific character- 
istics, of the environmental effects (decoherence) which 
can modify it [7]. In most cases, the influence of deco- 
herence during the evolution of a quantum walker trans- 
forms an originally pure state into a mixed state, low- 
ering the uncertainty about the location of the walker 
as it propagates. In the limiting case, when all cross- 
interference terms between different lattice sites are com- 
pletely erased, the state of pure diffusive classical propa- 
gation is reached [6]. 

QRWs have been theoretically explored for the case of 
one-dimensional lattices [5|, [8], and experimentally im- 
plemented by means of different physical platforms, such 
as photon-based systems (9l4T3|, optical lattices [l5| and 
waveguide arrays [16]. Also, QRWs have been imple- 
mented using trapped ions [17] and nuclear magnetic res- 
onance systems fl8l |. 

Although the implementation of one-dimensional 
QRWs has showed to be useful when describing several 
quantum information systems, there is a great interest 
in expanding the concept to multidimensional lattices. 



Along these lines, two-dimensional QRWs provide a pow- 
erful tool for modeling complex quantum information 
and energy transport systems fl9l [2p| . Notwithstand- 
ing, their realization represents a challenge because of 
the need of a four- level coin operation |2lH23l|. One way 
to overcome this drawback is to make use of different 
degrees of freedom of photons, such as polarization and 
orbital angular momentum, as it has been shown in |20(. 
Another approach is to mimic the two-dimensional QRWs 
evolution by performing two subsequent one-dimensional 
QRWs (HI!]!. 

Here, we make use of the latter approach to put for- 
ward an experimental setup for the realization of two- 
dimensional QRWs. We include the environmental effects 
(decoherence) as pure dephasing by means of the intro- 
duction of random phase patterns, generated by a spatial 
light modulator (SLM), which can be different from site 
to site (spatial disorder). By controlling the degree of 
decoherence, we study the transition from the quantum 
ballistic spreading to the diffusive classical walk. Also, 
by adding static disorder, we show the possibility of ob- 
serving Anderson localization [26] . 

Importantly, our proposal provides a versatile, highly- 
scalable experimental setup, which may be used as a tool 
for understanding quantum processes whose underlying 
physics can be somehow traced to the concept of ran- 
dom walks, such as energy transport in photosynthetic 
light- harvesting complexes [27|,l28| and material band gap 
structures [29] . 

The structure of the article is as follows. In Sec. II, we 
introduce the QRW model system, including dephasing, 
considered here. The proposed experimental setup is de- 
scribed in Sec. III. Numerical results are discussed in Sec. 
IV. Finally, we summarize our results in the conclusion. 



II. A TWO-DIMENSIONAL QUANTUM 
RANDOM WALK WITH DEPHASING 

A typical discrete quantum random walk comprises 
two operations: a coin-tossing operation and a shift op- 
eration. Here, the coin-tossing operation is performed 
in the Hilbert space % p spanned by vectors {\H), \ V}} : 
corresponding to the photon polarization. The random 
walk is performed in the Hilbert space Hx ® Hy, cor- 
responding to the position of the photon in the trans- 
verse plane, spanned by vectors {\i,j}} integers), 
which indicate sites (i, j) in the transverse plane (z, j = 
... — 2,-1,0,1,2...). The global quantum system thus 
evolve in the Hilbert space 

H = Hx®n Y ®n p . (1) 

The state of the system is described by the density 
matrix , which is transformed to a new density matrix 
each step n via the map 

^(n+l) = pW^^^H^t^^t^pWt, ( 2 ) 

H denotes the Hadamard operator 

*-£(!-.)• 

which acts on the polarization degree of freedom. The 
operators Sx and 5y, which describe the walker's shift 
in the transverse dimensions x and y, independently, read 
as 

S x = ^\i-lJ,H)(iJ,H\ + \i + lJ,V)(i,j,V\, (4) 
ho 

and 

Sy = M - l,H)(i,j,H\ + \i,j + l,V)(i,j,V\. (5) 

The coupling of the quantum walker with the environ- 
ment is described by pure dephasing [31] . The form of 
the unitary dephasing operator considered here can be 
written as 

& n) =Y,e-i+*?*'\i,j)(i,j\, (6) 

ij 

(n) 

where (p\- is a random phase matrix and <j z is the Pauli 

operator. Inspection of Eq. (|6j) shows that (f>ff repre- 
sents a newly introduced phase difference between the 
horizontal and vertical polarizations at each site. Con- 
cerning this, we will consider three physically relevant 
scenarios, that can be easily implemented in the set-up 
proposed here (see Sec. III). In the general case, the 

phase differences <p[^ are independent random variables, 
but with the same probability distribution. Moreover, 
the ensemble of phase differences <jffi can change from 
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FIG. 1. (Color online) General scheme for the implementa- 
tion of a two-dimensional random walk with decoherence. AT: 
attenuator; BS: beam splitter; HWP1 and HWP2: half-wave 
plates that act as coins in the random walk; HBD and VBD: 
horizontal and vertical beam displacers; M1,M2 and M3: mir- 
rors; SLM: spatial light modulator; CCD: spatial light sensor 
with single-photon sensitivity. 

step n to step n + 1. In the following, we will refer to this 
case as a QRW influenced by dynamical spatial disorder. 

The easiest probability distribution that we can con- 
sider is an uniform probability distribution. If phases can 
be chosen arbitrarily between the extreme values — £ and 
C, there is a constant probability l/(2£) to obtain any 
phase in this interval. £ = ir is the maximal phase which 
we can have between the two orthogonal polarizations. 
( = corresponds to the absence of any spatial disorder. 
If phases do not change during propagation, even though 
they might differ from site to site, i.e. — (j>^ +1 \ 

then we have static spatial disorder. Finally, if all phase 
differences are the same for all sites, but they can still 
change from one step to the following, we have dynami- 
cal dephasing without spatial disorder. 

The probability of detecting a photon in the site 

is 

P {n) (i,j) = (i,j\Tr p \pW]\i,j), (7) 

where the density matrix that describes the whole sys- 
tem is traced out over the polarization degree of freedom 
(Tr P ). 

The spreading of the uncertainty of photon location is 
characterized by dependence of the variance on the step 
index n 

V^=Y J P {n) {i,j)\n j -n\ 2 . (8) 

where r^- = (i,j) represents the lattice site with in- 
dexes and \i is the mean position, i.e., \i = 
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FIG. 2. (Color online) Detail of the SLM part of the setup 
allowing to generate dynamical spatial disorder. M3: mirror; 
EOD1 and EOD2: electro-optic deflectors; FL: Fourier lens; 
SLM: spatial light modulator; BS: beam splitter; IS1 and IS2: 
imaging systems 



the system. First, the polarization state of the pho- 
ton is changed by the half- wave plate (HWP1) to 
^|0, 0)0[(l + i)\H) + (1 - i)\V)], i.e., the Hadamard (H) 
operation is applied. After this transformation, the 
photon is displaced by the horizontal beam displacer 
(HBD) along the x axis according to its polarization, 
as described by the shift operator Sx- The photon is 
now in the state \ [(1 + i)\ - 1, 0, H) + (1 - i)|l, 0, V)] . 
A second half- wave plate (HWP2) implements a new 
Hadamard transformation, which transforms the quan- 
tum state to 4={(1 + i) (| - 1, 0, H) + | - 1, 0, V)) + (1 - 

i) (|1, 0, H) — |1, 0, V))}. The vertical beam displacer 
(VBD) shifts the position of the photon along the y axis. 
After this, the quantum state of the photon reads 



III. EXPERIMENTAL SETUP 

The main building block of the QRW setup is the 
multi-path Mach-Zehnder-like configuration shown in 
Fig. 1. It allows to make several runs of the QRW 
without the necessity of using a large amount of opti- 
cal components, as it is the case, for instance, of the 
experiment described in [Toj]. A similar scheme, based on 
a sin gle- path Mach-Zehnder-like configuration has been 
used PJ], LL4|]. However in these cases, the walker moves 
in time, whereas in our proposal the walker moves in 
the two-dimensional transverse plane, offering a way to 
simplify the experimental implementation of the two- 
dimensional QRW. 

As source of photons one can use a highly attenuated 
short coherent pulse, prepared by the combination of a 
photon source and attenuator (AT), generating the initial 
state 

\*M) = \0,0)®-L(\H)+i\V)), (9) 

where (0, 0) is the central site. The duration of the pulse 
has to be sufficiently smaller than the time-of-fly through 
the setup in one cycle. The transverse size of the Gaus- 
sian beam profile of the pulse has to be carefully cho- 
sen, so that two adjacent sites are not overlapping in the 
space due to diffraction. For instance, by making use of a 
Gaussian beam of 2 mm beam waist, corresponding to a 
Rayleigh range of 23.6 m (for a wavelength A = 532 nm), 
along with typical-sized optical components, we could in 
principle perform a QRW of approximately more than 
20 steps. The number of steps can be further improved 
by applying smaller beams together with a re-focusing 
system placed along the walker's path. Alternatively, a 
spontaneous parametric down-conversion source can be 
used, provided each down- converted photon is generated 
in a pure state. Then, the photon is transmitted via the 
beam splitter (BS) to the system. 

To get a clearer picture of the working of the quan- 
tum random walk, let us consider in detail the quan- 
tum state of the photon in its first passage through 



l*W) = -L[(i + i)(|-i,-i,ff) + |-i,i,F)) 

+ (l-i)(\l,-l,H)-\l,l,V))}. (10) 

The spatial light modulator (SLM) is used to introduce 
random phases, given by the phase-matrix 4^ • After 
the photon passes through the SLM, the whole cycle is 
repeated with high probability. If not, the photon escapes 
through the other output port of the beam-splitter, which 
directs it towards a highly sensitive (single-photon sen- 
sitivity) intensified CCD camera with integrated photo- 
multipliers, which allows to spatially resolve a weak sig- 
nal with a high efficiency (~90%) [32] . Moreover, losses 
in the setup can be compensated by an increase of the 
amplitude of pulse. The size of whole array of beams can 
be reduced by an auxiliary imaging system in order to fit 
on the limited size of sensitive area of CCD. 



A. Fast exchange of phase matrices 

A typical SLM has a response time in the order of 
tens of ms, which means that it is too slow for a fast 
phase-mask exchange. For this reason, the transmission 
SLM shown in Fig.l has to be supplemented by addi- 
tional components as it is shown in Fig. 2, which allows 
us to effectively generate the dynamical spatial disorder 
and dynamical dephasing without spatial disorder for a 
limited amount of steps. The time of exchange of phase 
matrices can be done now in tens of ps (33Ll33|, which 
is three-orders faster than it is indeed necessary in our 
experimental proposal, with a typical time of flight of 
the pulse in the order of ns. The scheme in Fig. 2 op- 
erates in the following manner. At the beginning, the 
size of the whole array of beams is reduced by an imag- 
ing system (IS1) to fit into the electro-optical deflector 
(EOD1). The EOD1, together with a Fourier lens (FL), 
serve to address different regions of the SLM along either 
the vertical or horizontal direction, which realize then the 
random phase-matrices (j>ff in all steps. Just changing 
the directions of deflection of the array, both types of 
dynamical disorders can be simulated. The EOD2 (and 
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another FL) then serve to return all beams back to the 
original direction of propagation, so that all beams re- 
main in the same position in the traverse plane at all 
time. 



IV. NUMERICAL RESULTS 



FIG. 4. (Color online) Spreading of the position of the pho- 
ton (V^) as a function of the number of steps for several 
values of £ and different types of dephasing. a) Dynamical 
spatial disorder, b) Dynamic dephasing without spatial dis- 
order. The results are obtained averaging over 500 different 
realizations of the matrix <j>ff . 



A. Quantum random walk 

Let us consider first the case when the SLM does not 
introduce any phase shift (cj)^ = for all This 
corresponds to the case £ = 0. Figure 2 shows the prob- 
ability distribution function p^ n \i : j) for (a) n = 10 and 
(b) n = 20 steps. In both cases, the distribution shows a 
symmetrical shape around the lines x = and y = 0, with 
four groups of peaks located along the x and y axes. The 
resulting symmetry comes from the specific initial quan- 
tum state chosen in Eq. (|9)) . When the number of steps is 
increased, the peaks move further away from the central 
site (0,0). The shapes obtained in Fig. 3 correspond to 
the probability distributions of a two-dimensional Grover 
walk [l9|,[25|]. The walker propagates with ballistic speed, 
characterized by a quadratic dependence of the variance 
with the step index, i.e., « n 2 . This case is 

shown in Figs. 4(a) and 4(b), corresponding to the case 
with C = 0. 



B. Quantum random walk affected by dephasing 

The dephasing effect introduced by the SLM allows 
to induce a transition from the quantum to the classical 
random walk via two mechanisms. First, as shown in 
Fig. 4(a), by means of dynamical spatial disorder. The 

(n) 

phase-matrix cj)\- ' shows independent and randomly cho- 
sen values for each site, and it is refreshed each step. 
The case £ = corresponds to the QRW with no dephas- 
ing. Increasing the amount of disorder, characterized by 
a corresponding increase of the parameter £, reduce the 
spreading of as can be seen in Fig. 4(a). In the 

limiting case, which is reached for £ = 7r, the observed 
dependence (~ n) of variance is a direct indication 
of the transition to the classical regime of random walks. 

The classical limit can also be reached by means of 
dynamical dephasing without spatial disorder, as it is 
shown in Fig. 4(b). The reduction of the uncertainty of 
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disorder, V» - 51.25 for ( = tt after 20 steps. On the 
contrary, for dynamical dephasing without spatial disor- 
der, we have ~ 98.45 under the same conditions. 
Indeed, the n-dependence of the typical deviation, char- 
acteristic of the classical regime, it is not yet reached 
after 20 steps, as is readily observed in Fig. 4(b). 



C. Anderson localization 

In the context of our discussion, the Anderson local- 
ization is the reduction of spreading of the uncertainty of 
the photon position [26|. We will demonstrate that this 
effect can also be observed in the set-up considered here. 
In [35], Anderson localization was observed in the trans- 
verse plane of a light beam passing through a crystal 
with random static fluctuations of the index of refrac- 
tion. Since the randomness in the index of refraction is 
affecting only the phase of the propagating beam, it is 
possible to imitate these phase-fluctuations with a SLM, 
under the consideration of static spatial disorder, since 
Anderson localization does not appear with dynamical 
spatial disorder. 

In Fig. 5(a) we present a typical profile of the phase- 
matrix 4^ij \ independent of n, which leads to beam lo- 
calization. Fig. 5(b) shows the corresponding probability 
distribution of the photon position for this specific phase 
profile. Notice that it contains a strong peak located in 
the middle of the lattice. The presence of Anderson lo- 
calization is confirmed in Fig. 5(c), where we show the 
averaged probability distribution function for ( = 7r, ex- 
hibiting an exponential suppression of probabilities for 
sites distant from the center. For the sake of clarity, we 
also plotted in Fig. 5(d) two cuts of the averaged proba- 
bility distribution along the X and Y axes, to highlight 
this feature. 
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FIG. 5. (Color online) Observation of the spatial Anderson 
localization, a) Example of a matrix <j>^ for £ = tt that leads 
to the Anderson localization, b) Corresponding probability 
distribution p^ n \i,j) after 20 steps, c) Averaged probability 
distribution over 500 realizations, d) Cuts of the data shown 
in Fig. 5(c) along the X and Y axes passing the site (0,0). 



the photon position is less dramatic than in the case with 
dynamical spatial disorder. For the dynamical spatial 



V. CONCLUSION 

We have put forward a new, highly scalable and eas- 
ily implemented experimental configuration to observe 
spatial two-dimensional random walks under a great va- 
riety of circumstances, by means of the implementation 
of two consecutive one-dimensional random walks. The 
proposal makes use of only a small amount of simple op- 
tics components, and allows us to simulate many different 
quantum systems and protocols based on the quantum 
random walk concept. Additionally, by carefully control- 
ling the amount and type of disorder present in the sys- 
tem, we have shown the effects of different environmental 
effects: dynamical spatial disorder, dynamical dephasing 
without spatial disorder and static spatial disorder. The 
last case drove us to the observation of Anderson localiza- 
tion. The control of environmental effects is paramount 
importance in nearly all quantum systems. In some cases, 
it is even crucial to understand the dynamics experimen- 
tally observe. For instance, in light harvesting complexes 
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[36(, the interplay between coherent evolution and noise 
is a critical ingredient, which allows to obtain a high ef- 
ficiency energy transport. 
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